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In  the  first  six  months  of  thin  ->roject  we  have  concentrated  our 
effort  in  the  following  areas: 

I.  Group  Development 

We  have  devoted  considerable  effort  to  recruiting  and  organizing 
our  research  group  for  this  project.  We  have  added  the  two  postdoctoral 
fellows  called  for  in  our  proposal.  They  are  Dr.  John  Shaner  and  Dr. 

Richard  Brandt.  We  are  extremely  fortunate  to  have  been  able  to  recruit 
these  people,  as  both  are  very  outstanding  young  physicists.  A  brief 
summary  of  their  qualifications  follows. 

Dr.  Shaner  received  his  B.S.  from  MIT  in  1964  and  his  Ph.D.  from 
the  University  of  California  at  Berkeley  in  1969  under  Professor  E.  L.  Hahn. 
He  held  a  postdoctoral  position  at  Berkeley  from  1969  until  August  1970 
when  he  joined  our  group. 

He  has  worked  in  ultrasonic  attenuation  in  addition  to  the  fields 
already  represented  by  his  publications: 

1)  "Mossbauer  Studies  of  Spin  Flip  in  Antiferromagnetic  Hematite", 
(with  N.  Blum,  A.  J.  Freeman  and  L.  Grodzins)  J.  Appl.  Phys.  36,  1136  (1965). 

2)  "Shock  Induced  Demagnetization  of  YIG",  (with  E.  B.  Royce) 

J.  Appl.  Phys.  32,  492  (1968). 

3)  "Zero  Field  Free  Precession  Measurements  of  the  Hyperfine  Inter- 
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action  in  Cs  Vapor",  (with  S.  A.  Miller  and  E.  L.  Hahn)  Bull.  Am.  Phys. 
Soc.  13,  31)0  (1968). 

}i )  "Cyclotron  Echoes  in  a  Weakly  Ionized  Cesium  Plasma",  (with 
E.  L.  Halm)  J.  Appl.  Phys.  4l,  8j8  (1970). 

He  came  enthusiastically  recommended  by  Professor  Hahn  and  many  other 
physicists  who  have  had  the  opportunity  to  observe  him.  Since  arriving  lie 
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has  thoroughly  lived  up  to  our  expectations. 

Dr.  Brandt  received  his  B.S.  from  Cal  Tech  in  1963  and  his  Ph.D. 
from  the  University  of  Illinois  in  1967  under  Professor  Fred  Brown.  From 
1967  to  the  present  he  has  been  at  Lincoln  Labs  working  primarily  on  far 
infrared  spectroscopy.  He  will  Join  us  early  in  1971* 

His  publications  include: 

1)  "Etching  on  High  Purity  Zinc",  (with  K.  H.  Adams  and  T.  Vreelaud,  Jr} 
J.  Appl.  Phys.  34,  587  (1963)- 

2)  "Dislocations  and  Etch  Figures  of  High  Purity  Zinc",  (with  K.  II. 

Adams  and  T.  Vreeland,  Jr.)  J.  Appl.  Phys.  ^4,  591  ( 19^3 ) • 

3)  "Infrared  Absorption  and  Vibronic  Structure  due  to  Localized 
Polaions  in  the  Silver  Halides",  (with  F.  C.  Brown)  in  Localized  Excitations 
in  Solids,  edited  by  R.  F.  Wallis  (Plenum  Press,  Inc.,  New  York,  1968). 

4)  "Polaron  Zeeman  Effect  in  AgBr",  (with  D.  M.  Larsen,  P.  P.  Crooker 
and  G.  B.  Wright)  Phys.  Rev.  Letters  2£,  240  (1969)* 

5)  "Induced  Infrared  Absorption  due  to  Bound  Charge  in  the  Silver 
Halides",  (with  F.  C.  Brown)  Phys.  Rev.  iBl,  124l  (1969). 

6)  "A  Beam  Condenser  for  Infrared  Spectrophotometers",  Appl.  Optics 
3,  315  (1969). 

7)  "Polaron  Zeeman  Effects  in  the  Silver  Halides",  (with  P.  M. 

Larson  and  D.  R.  Cohn)  Proceedings  of  the  10th  International  Confeivneo 
on  Physics  of  Semiconductors,  Cambridge,  August  1970.  to  be  published. 

Dr.  Brandi,  was  enthusiastically  recommended  by  many  physicists 
■>ic lading  Piwfessors  Muiurr and  Brown  at  Illinois.  However,  the  strongest 
rveomiiiradul.ioii  of  al  l  i that  he  turned  down  a  posit  ion  as  Assistant 
Professor  al  the  University  of  Illinois  to  Join  our  group. 

We  mil  ieipafe  Hr.  Hmndl.  In* lug  a  givul.  asset  l.o  the  group,  pai  l. irn |.ar l,\ 


since  he  is  an  established  expert  in  far  infiared  spectroscopy,  and  has 
had  considerable  experience  in  research  on  wide-gap  semiconductors. 

In  addition  to  Brandt  and  Shaner  whom  we  have  brought  into  the 
Department  we  have  enlisted  the  help  of  some  people  already  on  the  staff. 
Professor  Frank  Harris  has  become  interested  in  an  empirical  band  structure 
calculation  on  rutile  and  has  donated  the  services  of  one  of  his  postdoc's 
to  the  calculation— under  his  supervision.  Professor  Harris  has  an 
international  reputation  as  a  theorist  and  we  are  very  pleased  to  have  his 
help.  Results  of  this  calculation  should  be  very  helpful  in  interpreting 
electronic  measurements. 

We  also  have  two  graduate  students — one  working  on  photoconductivity 
and  one  on  crystal  growth. 

II.  Crystal  Growth 

A  good  part  of  the  crystal  growth  program  during  the  first  six  months 
under  this  grant  has  involved  a  comprehe naive  literature  study  of  potential 
crystal  growth  techniques  for  rutile  and  the  physical  chemistry  of  certain 
titanium  compounds  which  might  be  used  as  starting  materials  for  the  prepa¬ 
ration  of  ultrapure  TiOg  in  later  stages.  Consequently  detailed  plans  for 
the  crystal  growth  and  purification  equipment  as  well  as  consecutive  goals 
in  the  procedures  were  worked  out.  The  crystal  growth  furnace,  a  novel 
development,  has  been  built  and  tested  already.  First  growth  runs  will 
be  performed  within  the  next  two  months. 

III  .Dielectric  Loss  Measurements 

Current  status  of  this  part  of  the  project  is  summarized  in  Appendix 
I,  which  is  a  preprint  of  a  paper  to  be  submitted  shortly. 

IV  Dielectric  Loss  Analysis 

This  aspect  of  the  program  has  turned  out  to  be  a  rather  formidable 
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problem  and  has  required  a  substantial  effort.  The  problem  is  best,  handled 
in  two  steps;  the  dielectric  loss  and  capacitance  spectrum  are  first 
"inverted"  to  obtain  a  resistivity  spectrum,  which  must  then  be  ordered 
in  such  a  way  as  to  satisfy  Poisson's  equation.  The  first  part  of  the 
procedure  lias  now  been  reduced  to  a  routine  computer  analysis  and  we  have 
made  c insiderable  progress  in  handling  the  second  step;  further  progress 
will  require  knowledge  of  electron  trapping  levels,  which  we  expect  to 
obtain  from  other  measurements.  The  details  of  the  analysis  of  dielectric 
loss  data  are  also  summarized  in  Appendix  I. 
y.  Review  Paper 

Probably  because  of  its  unique  properties,  TiO^  has  attracted  a 
great  deal  of  interest  and  research  is  being  carried  on  by  investigators 
in  many  different  disciplines.  Because  of  this,  research  results  are 
published  in  a  wide  variety  of  Journals,  many  of  which  are  not  readily 
accessible.  The  extent  of  the  research  activity  in  this  field  is  quite 
startling;  in  the  past  three  years,  more  than  250  technical  papers  concern¬ 
ing  the  physical  properties  of  rutile  and  closely  related  compounds  have 
been  published.  We  felt  it  was  absolutely  essential  to  thoroughly  digest 
all  of  this  material  before  proceeding  very  far  on  most  of  our  investi¬ 
gations,  particularly  since  some  of  our  group  had  not  previously  worked 
on  this  material.  Any  group  working  in  this  area,  of  course,  is  (or 
should  be)  confronted  by  the  same  problem.  The  last  comprehensive  review 
paper  in  this  area  was  published  in  1959?  &  really  thorough  and  critical 
review  is  badly  needed.  We  felt  wc  could  perform  a  substantial  service 
to  the  scientific  community  by  summarizing  the  results  of  our  extensive 
literature  search  and  evaluation  in  the  form  of  a  review  article.  We 
have  already  received  the  blessing  of  the  editors  of  Reviews  of  Modern 
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Physics  for  this  project  and  are  well  on  our  way  to  completing  this  work. 

We  have  on  file  virtually  all  of  the  significant  papers  published  in  this 
area  in  the  past  ten  years,  have  cataloged  and  cross-referenced  these  (see 
attached  outline)  and  have  already  completed  evaluation  of  many  of  them. 

We  have  begun  writing  the  sections  concerning  electrical,  optical,  and 
mechanical  properties,  and  the  phonon  spectrum. 

This  project,  in  addition  to  serving  a  very  useful  purpose  for  other 
workers  in  the  field,  has  been  extremely  helpful  to  us  directly,  in  that 
it  has  brought  out  the  importance  of  several  areas  of  investigation  and 
has  suggested  several  critical  experiments  which  should  help  to  resolve 
some  of  the  apparent  contradictions  in  the  literature.  In  particular, 
optical  measurements  on  conduction  electrons  at  low  te rape rat  ores,  far 
infrared  spectroscopy  for  study  of  defects  and  ultrasonic  attenuation 
measurements  look  very  promising  (these  are  discussed  in  detail  in  Sec.  II). 

In  addition,  we  have  set  up  the  experimental  equipment  to  do 
photoconductivity  measurements  and  have  made  several  runs  with  it.  There 
do  not  appear  to  be  any  serious  experimental  problems,  but  the  results 
are  very  complex  and  we  have  not  yet  achieved  sufficient  understanding 
to  merit  a  detailed  discussion. 

We  have  nearly  completed  the  experimental  set  up  for  photoeniission 
studies  but  have  not  yet  taken  any  actual  data. 

Our  work  on  devices  thus  far  has  consisted  of  formulating  some 
specific  configurations  we  would  like  to  investigate.  Potential  problems 
associated  with  fabrication  of  these  devices  have  been  discussed  with 
IV.  Kobe r l  Huber  of  General  Instruments,  ami  arrangements  for  production 
of  tee  devices  by  General  Instruments  iiave  been  completed. 


Working  Outline  for  the  Rutile  Review  Paper 


I.  Crystal  Growth  and  Preparation 

a.  Structure  and  Stoichiometry 

/ 

b .  Growth 

c.  Purification  / 

d.  Sample  Preparation 

II.  Electrical  Properties 

a.  Conductivity 

b.  Hall  Effect  Measurements 

c.  Photoconductivity 

d.  Piezo- resistance 

e.  Dielectric  Loss 

f.  Ultrasonic  Attenuation 

III.  Thermal  Properties 

a.  Thermal  Conductivity 

b.  Thermoelectric  Effects 

c.  Specific  Heat 

d.  Anharmonic  Effects 

IV.  Optical  Properties 

a.  Optical  Absorption  and  Reflection 

1.  Ultraviolet 

2.  Visible 
3*  Infrared 

b.  Modulation  Spectroscopy 
o.  Dielectric  Function 

n.  Hainan  Scattering 

o.  Pliotoemission 

f.  Thermal  Emission  and  Related  Effects 


V. 


Magnetic  Properties 


a.  Susceptibility 

b.  Nuclear  Magnetic  Resonance  and  Mttssbauer  Effect 

c.  Electron  Paramagnetic  Resonance 

d.  Maser  Action 

VI.  Electronic  Structure 

a.  Intrinsic  Band  Structure 

b.  Defect  and  Impurity  Levels 

VII.  Phonon  Spectrum 

VIII.  Intrinsic  Defects  (Ti  or  0) 

IX.  Impurities  (Defects  other  than  Ti  or  0) 

X.  Mechanical  Properties 

a.  Internal  Friction 

b.  Elasticity  Experiments 

XI.  Radiation  Damage 

XII.  High  Pressure  Properties 
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APPENDIX  I 


Dielectric  Loss  Mechanism  in  Rutile  (TIOg)* 

0.  W.  Johnson  and  John  W.  DeFord 
Department  of  Physics,  University  of  Utah 
Salt  Lake  City,  Utah  84112 

Sverre  ttyhra 

Department  of  Physics,  University  of  Queensland 
Brisbane,  Australia 

\ 

A  new  mechanism  to  account  for  the  anomalously 
large  dielectric  loss  and  apparent  dielectric 
constant  frequently  observed  in  TiOg  is  presented, 
which  appears  to  agree  well  with  experiment.  The 
mechanism  involves  field- induced  donor  migration, 
resulting  in  a  large  decrease  in  electronic  resis¬ 
tivity  in  part  of  the  crystal.  Experimental 
results  on  a  variety  of  crystals  are  discussed, 
and  the  influence  of  electrode  materials,  surface 
condition,  heat  treatment  and  ambient  is  analyzed. 
A  method  for  numerically  calculating  the  resis¬ 
tivity  profile  of  such  a  crystal  from  dielectric 
loss  data  is  described.  Results  indicative  of 
room  temperature  p-type  conductivity,  and  a 
voltage-induced  conduction  process  are  presented. 


*Work  supported  by  AREA,  Grant  No.  DAHC  15  G  13- 


I.  Introduction 
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Previous  dielectric  loss  studies  in  rutile  have  yielded  results 

which  have  been  interpreted  to  imply  anomalously  large  frequency  dependent 

4 

dielectric  constants  ranging  as  high  as  10  ,  both  for  oxidized  crystals  and 

slightly  reduced  crystals,  in  which  two  terminal  dc  resistance  measurements 

13  8  12 

indicate  a  resistivity  ~  10  ohm-cm  for  oxidized,  and  from  10  to  10 

ohm-cm  for  reduced  crystals.  In  contrast,  microwave  measurements1*1'0  have 

-12 

consistently  yielded  e  «  173  eQ  (gq  =  8.85  x  10  in  MKB  units)  in  the 
c-axis  direction.  The  anomalously  large  apparent  dielectric  constants  have 
generally  been  accompanied  by  dielectric  loss  peaks  (D  =0.5  or  more,  even 

3 

in  heavily  oxidized  crystals),  typically  in  the  frequency  range  from  10  to 

4 

10  Hz.  A  layer  of  material  of  resistivity  p  will  exhibit  large  electronic 

dissipation  only  at  frequencies  less  than  ~  l/(2TTpe)*^  Assuming  e  «  170  eQ 
13 

and  p  w  10  fi  cm,  significant  electronic  dissipation  should  be  observed 

.3 

only  below  ~  10  Hz.  Thus,  if  the  loss  peaks  are  electronic  in  origin 

large  dissipation  in  the  kilocycle  range  would  require  either  material 
6  7 

with  p  of  10  or  10'  fi  cm,  or  a  very  large  and  strongly  frequency  depen¬ 
dent  dielectric  constant.  Mechanisms  to  account  for  these  results  have 
generally  involved  a  "depletion  layer"  (i.e.,  a  layer  in  which  the  "normal" 
concentration  of  conduction  electrons  has  been  reduced),  induced  by  the 
metal  electrodes  in  contact  with  the  crystal,  or  an  enhancement  of  the 
dielectric  constant  due  to  the  presence  of  conduction  electrons.  As  we 
shall  show,  the  first  mechanism--by  itself— will  not  produce  the  observed 
dielectric  loss  peaks.  The  second  mechanism  is  readily  ruled  out  by 
calculating  the  contribution  of  conduction  electrons  to  the  real  part 
of  the  dielectric  constant. 11  It  is  many  orders  of  magnitude  too  small 
to  account  for  the  observed  effects.  Motion  of  any  other  charged  species 


present  could  contribute  to  dielectric  loss,  but  none  of  the  known 
defects  or  impurities  in  rutile  diffuse  rapidly  enough  at  room  tempera¬ 
ture  to  contribute  significantly.  Other  loss  mechanisms,  such  as  dipolar 
reorientation,  motion  of  charge  dislocations,  loss  due  to  conductive 
inclusions,  etc.,  also  appear  to  be  ruled  out  by  the  known  properties 
of  the  crystal. 


II.  Analysis  of  Electronic  Dielectric  Loss 
In  the  absence  of  any  other  plausible  loss  mechanism,  we  have 
attempted  to  account  for  the  observed  behavior  on  the  basis  of  changes 
in  bulk  resistivity  of  the  crystal,  induced  by  contact  with  the  electrode 
material.  We  assume  tnat  none  of  the  bulk  parameters  of  the  crystal, 
except  for  the  resistivity,  are  changed  by  the  presence  of  the  electrodes, 
and  that  the  thickness  of  the  crystal  is  small  compared  to  its  lateral 
dimensions.  The  required  resistivity  spectrum  can  be  derived  from  the 
experimental  measurements  in  a  fairly  straightforward  way.  By  "resistivity 
spectrum"  we  mean  the  total  thickness  of  layers  having  resistivity 
The  analysis  gives  no  information  as  to  the  "resistivity  profile",  i.e., 
the  relative  position  of  the  layers.  The  resistivity  profile  can  in 
principle  be  obtained  by  requiring  that  the  layers  be  arranged  to  satisfy 
Poisson's  equation  with  appropriate  boundary  conditions. 

We  can  derive  the  resistivity  spectrum  as  follows:  the  complex 


impedance  z  of  a  crystal  of  thickness  l  (assuming  planar  geometry)  is 

2-  g&riwj dx 

0J  L  p2(x)(j)  e  +  1  J 


=  1  d6 
oJ  p  (eke  +  1 


(1) 
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where  <d  is  the  angular  frequency.  The  experimentally  measured  parameters 
are  the  series  equivalent  capacitance  C  and  the  dissipation  D,  given  by 


C  = 


djlm(z) 


and 


Im(z) 


Thus,  z(ou)  may  be  determined  from  the  experimental  parameters.  Our 

problem  is  to  solve  the  integral  equation  (l)  for  p.  There  are  several 

possible  approaches  to  the  problem,  but  all  involve  approximating  the 

integral  by  a  sum:  2 

pv  -  ipv  m.e 

?<«,  )  „  i  ?  Js - S-J, 


•J' 


(p^j* 


p7T9k’^i Vk  (2) 


where  is  the  average  resistivity  of  the  k  layer,  of  thickness  -20^. 

In  this  form  it  is  obvious  that  the  ordering  of  the  layers  is 
immaterial  find  hence  that  only  the  resistivity  spectrum  can  be  determined 
in  this  fashion.  One  possible  approach  at  this  point  is  simply  to  assume 
a  set  of  p's.  This  determines  the  a^  in  Eq.  (2).  If  n  p's  are  chosen 
we  can  then  use  (n-l)  data  points  to  get  (n-l)  equations  in  the  n  unknowns, 
0^.  By  then  imposing  the  condition 


we  get  a  system  of  n  equations  in  n  unknowns  which  can  be  solved  exactly. 
Uhf ort unate ly  the  exact  solutions  obtained  in  this  way  are  generally 
non-physical  (some  of  the  0fc  turn  out  to  be  negative),  due  to  the  approx¬ 
imations  inherent  in  the  technique  and  small  inaccuracies  in  the 
experimental  data.  Thus,  a  somewhat  more  sophisticated  treatment  is 
required. 

The  first  improvement  is  to  use  a  form  of  least  squares  fit  rather 
than  attempting  an  exact  solution  to  an  inexact  problem.  Instead  of 
z(uj)  it  is  convenient  to  use  the  dissipation  D(uj)  (either  C( ou)  or  D( ou) 


may  be  used  since  they  are  related  by  Kramers -Kronig  equations).  Approx¬ 
imating  the  integrals  by  sums  we  have 

Pk9k 


n 


£  (p^ise)  +  1 


D(oj)  = 


k=l 


n 

kil 


2, 

pk  ekeco 


(3) 


(pkuje)£ 


+  1 


where  we  nov  restrict  the  fractional  thicknesses,  0^,  to  non- negative 

values;  n  is  the  assumed  number  of  layers.  A  reasonably  good  (~  2$) 

fit  to  the  experimental  D(uj)  can  be  found  by  trying  various  p  and  0  choices, 

plotting  out  the  resulting  theoretical  curves  (calculated  from  Eq.  (3) ), 

and  adjusting  the  pk  and  0^  to  get  the  best  fit.  This  method  works 

reasonably  well,  because  a  given  resistivity  layer  affects  D  only  over 

g 

a  fairly  narrow  frequency  range  and  hence  the  corrections  in  p  can  be 
made  more  or  less  systematically.  Unfortunately,  it  is  tedious,  time 
consuming,  and  not  suitable  for  routine  analysis.  Hence  we  have  developed 
a  more  systematic  approach  using  the  1108  computer.  This  approach  is 
described  in  Appendix  A.  It  leads  to  a  resistivity  spectrum  which  gives 
the  best  least  squares  fit  between  calculated  and  experimental  dissipation 
curves. 

Assuming  the  resistivity  spectrum  is  now  known,  the  corresponding 

conduction  electron  density  is  determined,  in  terms  of  the  electron 

mobility  |a.  Knowledge  of  the  density  of  states  in  the  conduction  band 

would  then  permit  calculation  of  the  position  of  the  bottom  of  the 

conduction  band,  E,  relative  to  the  Fermi  level,  E_.  We  now  formulate 

F 

the  problem  of  ordering  the  layers  to  obtain  the  resistivity  profile. 


Poisson's  equation  may  be  written 

=  -  f  =  -  f  ["+  -  -  nJ  w 

where  ♦  is  the  electrostatic  potential,  p  is  the  net  charge  density, 
e  =  1.6x10  ^  coul.,  n,j,  and  nQ  are  the  densities  of  electrons  in  traps 
and  the  conduction  band,  respectively,  and  eN+  is  the  net  positive  charge 
density,  due  to  impurities,  defects,  etc.  We  assume  that  before  applica¬ 
tion  of  electrodes,  ♦  is  constant  through  the  crystal,  so  that  N+  =  +  nQ. 

We  assume  throughout  that  is  well  above  the  valence  band  so  that  hole 
density  can  be  ignored  (although  there  is  some  evidence  that  this  assumption 
may  not  be  valid  in  certain  cases,  as  discussed  in  the  next  section). 

will  consist  of  various  substitutional  and  interstitial  species,  some 
of  which  apparently  are  mobile  at  temperatures  of  ~  300°C;  it  is  convenient 
to  write 

"+  -  "f  +  f  ZI  «I 

where  Nj,  includes  all  immobile  ions,  and  eZ^  is  the  charge  of  the  Ith 
mobile  ion  species,  with  density  N^. 

In  the  absence  of  space  charge,  Ec  is  just  the  binding  energy  of 
a  conduction  electron  to  the  crystal,  and  is  fixed  relative  to  the  vacuum 
level.  Space  charge  layers  in  the  crystal  result  in  a  shift  of  the 
potential  energy  of  an  electron  in  the  conduction  band  in  the  interior  of 
the  crystal;  thus,  in  general  Ec  will  be  a  function  of  position,  and  it 
is  clear  that 

Ec(x)  =  -  e4»(x)  +  Y 

where  y  is  a  constant  which  depends  on  our  choice  of  zero  of  potential. 
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It  is  also  clear  that  Ec  at  the  surface  of  the  crystal  is  constant  relative 

to  the  vacuum  level  in  the  absence  of  space  charge  external  to  the  crystal. 

We  may  write  expressions  for  nQ  and  n^,  in  terms  of  Ec  and  the 

Fermi  energy  E_: 

r  oo 

„  =  , _ 

c  Er  e^E  EF^  +  1  eE  Ec'EF  +  1 
c 

"T  =  i  efKVV  +  ! 

where  D(E)  is  the  density  of  states  in  the  conduction  band  (for  which 
we  use  the  usual  approximation  of  an  effective  density  of  states  nQ 
located  at  Ec),  p  =  l/luT,  and  n^  is  the  density  of  traps  with  trapping 
energy  E^. 

Thus  Eq.  4  becomes 


Next  consider  the  situation  with  electrodes  applied.  At  equilibrium, 
the  Fermi  level  will  be  constant  throughout  the  system.  If  mobile  ions 
are  present,  a  similar  condition  must  be  satisfied  for  the  ions  as  well — 
that  is,  the  chemical  potential  of  each  mobile  species  must  be  constant 
throughout  the  system.  Before  the  electrodes  are  applied  we  have  the 
situation  shown  in  Fig.  1A. 

When  the  electrodes  are  applied,  charge  will  flow  until  E_  is 
constant  throughout  the  system.  Since  the  density  of  states  is  much 
higher  in  the  metal  than  in  the  rutile,  the  Fermi  level  in  the  metal 
will  not  change  appreciably.  Furthermore,  the  resulting  space  charge 

.7 

in  the  metal  will  be  confined  to  a  very  thin  layer  (<  10  cm)  adjacent 
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to  the  rutile;  hence,  Ec  at  the  surface  of  the  rutile  is  essentially 
unchanged  by  the  presence  of  the  electrodes.  There  are  three  possibilities: 
1)  The  Fermi  level  is  initially  higher  in  the  metal  than  in  the  rutile, 
i.e. ,  the  work  function  of  the  rutile  is  >  the  work  function  of  the  metal. 
(We  refer  throughout  to  the  thermodynamic  work  function,  the  difference 
between  the  vacuum  level  and  Ec* )  2)  The  Fermi  level  is  initially  lower 

in  the  metal  than  in  the  rutile.  3)  The  Fermi  level  is  initially  the 
same  in  both  materials.  In  case  l),  electrons  will  flow  from  the  metal 
to  the  rutile  raising  Ep,  in  the  rutile,  resulting  in  a  net  positive  charge 
in  the  metal  at  the  interface  and  a  net  negative  charge  in  the  rutile.  Due 
to  the  planar  geometry  we  can  use  Gauss's  Law  to  determine  the  fields  and 
energies  involved.  We  can,  without  significant  error,  treat  the  metal 
space  charge  as  a  charge  layer  a+  at  x  =  o,  the  rutile-metal  interface. 

Then  the  change  in  electrostatic  potential,  A<1>,  at  a  depth  x  in  the  rutile 
is  given  by 

X  X  x' 

AO  =  E(x,)dx'  *=  -7  ^  [o+  ^x'^dx'^dx’ 

O  -  O-O' 

where  p”  is  the  net  charge  density  in  the  rutile,  and  E  is  the  resulting 

electric  field.  Thus,  A*  -*  o  as  x  "*  o,  so  that  the  energy  levels  at  the 

surface  of  the  rutile  are  unchanged  by  application  of  electrodes,  as 

x 

previously  noted.  Note  also  that  ^  {r(x)  dx  -*  o+  as  x  -*  ®,  in  order  to 

o 

o 

maintain  overall  charge  neutrality.  The  electron  energy  levels  in  the 

interior  are  shifted  by  the  electrode-induced  space  charge  by 

x  x' 

AEc  =  (ea*)  =  |  ^  [ct+  -  \  p(x")dx"l  dx’. 

O-  cr  J 

Hence,  the  resulting  configuration  is  as  indicated  in  Fib.  IB.  The 
thickness  of  the  space  charge  layer  can  be  calculated  numerically  from 


the  initial  difference  in  Fermi  levels,  if  the  trap  density  is  known 
and  there  is  no  positive  ion  migration.  For  physically  reasonable  dif¬ 
ferences  in  work  function,  the  calculated  thickness  of  the  layer  is 
always  <  Ip,  even  in  the  absence  of  electron  traps  (traps  reduce  the 
thickness).  Hence,  the  effect  of  the  electrodes,  in  the  absence  of 
ion  motion,  is  to  induce  a  very  thin  (typical  sample  thickness  is  ~  500p) 
conducting  layer  near  the  surface,  leaving  the  bulk  of  the  crystal 
unchanged. 

Curve  A  in  Fig.  2  presents  the  results  obtained  from  a  numerical 

integration  of  Eq.  (5),  assuming  Eq  -  Ep  before  electrodes  were  applied 

was  —  1.5  eV  at  room  temperature,  that  conduction  band  density  of  states 

20,  3 

could  be  represented  by  an  effective  density,  nQ  of  10  /cm  at 
E  ,  and  that  the  difference  in  work  functions,  (W  -W  )  >  1.35  eV.  We 
further  assumed  that  there  were  no  electron  trapping  levels  above  E_; 
this  is  probably  not  the  case,  but  failure  of  this  assumption  would  not 
result  in  major  changes  in  the  calculated  parameters.  Our  results  are 
not  particularly  sensitive  to  our  choice  of  nQ,  either,  except  for  the 
calculated  values  of  donor  density  discussed  below,  which  must  be  con¬ 
sidered  rough  estimates.  For  the  assumed  parameters,  D  is  negligible 
and  C  constant  over  our  frequency  range,  as  indicated  by  Curves  A  of 
Fig.  3* 

The  question  of  possible  ion  motion  must  now  be  considered.  It 
is  readily  shown  that  the  effect  of  a  field  E  on  a  charged  impurity  is 
to  produce  a  drift  velocity 


ZeE  _ 

Va“kf-D 

where  D  is  the  thermal  diffusion  coefficient,  and  Ze  is  the  ion  charge. 
While  precise  values  of  D  are  not  available  for  any  of  the  donor  irn,  'ities 
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(except  Li,  which  is  probably  not  present  in  appreciable  concentrations) 
reasonable  estimates  indicate  that  many  interstitial  ions  should  move 
significant  distances  at  temperatures  of  ~  300°C,  in  the  presence  of  a 
field,  even  though  significant  thermal  diffusion  is  not  observed  except 
at  substantially  higher  temperatures.  In  particular,  H+  is  certainly 

mobile  in  this  range  (as  discussed  below)  and  other  species  with  ionic 

♦  4+  13  l4 

radius  <  0-70 A,  including  Ti  ,  may  move  to  a  lesser  extent.  ’ 

Thus,  we  assume  that  there  are  positive  ions  present  which  are 
mobile  in  the  presence  of  an  electric  field,  at  least  at  temperatures 
of  the  order  of  200  -  300°C,  as  indicated  by  experimental  results  described 
in  Sec  III.  In  the  electron  space  charge  region  these  ions  are  subjected 
to  a  large  electric  field  (>  10^v/cm),  driving  them  away  from  the  electrodes, 
due  to  the  combined  effect  of  the  net  positive  charge  at  the  interface 
and  the  distributed  rutile  space  charge  layer.  This  will  tend  to  remove 
positive  ions  from  a  layer  beneath  the  electrodes  until,  at  equilibrium, 
the  ion  concentration  gradient  just  balances  the  electric  field. 

It  is  straightforward  to  show  that  at  equilibrium,  the  concentration 
of  ions  is  given  by: 


where  Zj.e  is  the  charge  of  the  ion,  Nq  and  are  the  values  at  x  =  0,  and 
♦  (x)  is  the  electrostatic  potential  at  x.  (After  ion  migration — the  change 
in  ion  distribution  alters  ♦  from  the  value  previously  calculated,  of  course.) 
It  follows  that  the  fractional  change  per  unit  distance  is  given  by: 

1  dN  ZeE 
N  dx  ”  kT 

where  E  =  -d$/dx.  Thus,  the  concentration  will  change  with  a  "characteristic 
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length"  given  by 


Ax 


kT 

ZeE 


A  space  charge  layer  containing  10^  electrons/ cm^  results  in  &  Mi  ■ 

As  shown  below,  a  typical  space  charge  layer  has  roughly  10"^  electrons/ cm^ 
hence,  the  concentration  will  change  over  a  distance  which  is  small  com¬ 
pared  to  the  typical  sample  thickness,  resulting  in  an  equilibrium 
positive  ion  distribution  which  is  essentially  a  step  function.  Note 
that  it  is  possible  for  this  depletion  layer  to  be  much  thicker  than  the 
original  space  charge  layer,  since  the  region  of  high  field  "migrates" 
into  the  crystal  as  the  positive  ions  are  displaced. 

It  is  readily  seen  that  the  central  region  of  the  crystal  is  in 
a  state  of  almost  exact  charge  neutrality,  since  any  significant  departure 
from  neutrality  results  in  further  ion  migration  and  hence  does  not 
represent  equilibrium.  An  electric  field  is  possible  only  in  the  narrow 
region  where  the  positive  ion  distribution  is  changing,  and  in  the  ion 
depletion  layer.  The  exact  configuration  may  be  obtained  from  Eq.  (5) 
with  appropriate  substitutions  for  the  various  change  density  terms, 
and  with  the  boundry  condition  on  Ec(0)  fixed  in  terms  of  the  relative 
work  functions  of  the  rutile  and  the  electrode.  A  subsidiary  condition 
on  the  total  number  of  mobile  ions  must  be  applied  to  obtain  a  unique 
solution  to  this  equation.  In  the  case  where  only  ions  already  in  the 
crystal  are  involved,  the  condition  is  simply  that  the  total  number  of 
ions  in  the  crystal  must  be  constant.  If  donor  impurities  may  diffuse 
in  from  the  electrodes  or  the  ambient,  the  final  concentration  is  fixed 
from  thermodynamic  considerations  and  represents  the  solubility  limit 
for  the  given  experimental  conditions  (this  point  is  discussed  in  detail 
in  the  next  section).  For  the  assumed  values  of  the  parameters  of  the 
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specimen  described  above,  the  result  of  numerical  integration  of  Eq.  (5) 
for  a  temperature  of  350°C  is  shown  in  Fig.  2,  Curve  B.  For  this  cal¬ 
culation,  an  initial  density  of  mobile  ions  of  4x10^/ cm^  was  assumed, 
and  it  was  assumed  that  no  donor  impurities  diffused  in  from  the  ambient. 
This  density  is  probably  too  low  but  essentially  indistinguishable  results 
obtain  with  a  higher  density  of  mobile  ions,  so  long  as  they  are  compen¬ 
sated  by  electron  trapping  states,  as  would  be  required  to  produce  the 
high  initial  resistivity.  For  the  assumed  initial  density,  equilibrium 
at  350°C  requires  essentially  complete  depletion  to  a  depth  of  ~  6^, 
then  very  rapid  increase  to  a  constant  value  of  ~  4 .  lxlO1^/ cm^  • 

The  behavior  indicated  in  Fig.  2  is  a  general  feature  of  all  of 
our  calculations;  it  turns  out  that  for  any  reasonable  choice  of  para¬ 
meters  and  for  thicknesses  of  >  10C|i,  there  is  always  a  broad  central 

/v 

region  of  essential  charge  neutrality  and  uniform  ion  density.  Even 

in  q 

small  net  charge  densities  (~10  /cnr)  quickly  cause  E  to  diverge 

beyond  physically  realizable  limits. 

When  the  specimen  is  quenched,  the  ion  distribution  will  presumably 

be  "frozen  in"  (we  assume,  as  indicated  by  experiment,  that  the  ions  are 

relatively  immobile  at  room  temperature),  whereas  the  electron  distribution 

changes  to  re-establish  a  constant  E_,.  However,  it  must  be  recognized 

£ 

that  the  electron  concentration  near  the  center  of  the  specimen  remains 
essentially  at  the  value  established  at  the  annealing  temperature,  in 
order  to  preserve  charge  neutrality,  as  discussed  above.  This  is  the 
crucial  point  in  our  analysis,  since  this  requires  that  the  quantity 
[(Ec  -  EyVkT]^^  be  the  same  before  and  after  quenching  to  room  temperature, 
resulting  in  a  substantial  lowering  of  Ec(l/2).  (Note  that  the  presence  of 
unfilled  electron  traps  would  reduce  the  magnitude  of  the  effect  but  would 


not  alter  the  essential  conclusions  of  the  analysis.)  Again  solving 
Eq.  (5)  for  our  assumed  parameters,  but  with  the  ion  distribution  calculated 
in  Curve  B  held  fixed,  results  in  Curve  C  of  Fig.  2.  This  result  is  rather 
surprising;  we  see  that  the  redistribution  of  the  positive  ions  has  resulted 
in  a  wide  region  of  relatively  low  resistivity  in  the  interior,  with  a 
thin  layer  of  high  resistivity  material  under  the  electrodes.  As  will  be 
seen  in  the  following  section,  this  is  precisely  the  type  of  resistivity 
profile  required  to  account  for  the  experimental  data.  The  calculated 
values  of  C  and  D  for  this  specimen  are  shown  as  Curves  B,  Fig.  3-  The 
upward  curvature  of  Ec(x),  moving  from  the  interior  toward  the  surface, 
results  from  a  very  slight  departure  from  charge  neutrality  (electron 
deficiency)  in  the  center;  this  is  the  only  conduction  band  configuration 
which  converges  to  the  correct  value  at  the  surface,  with  the  calculated 
positive  charge  distribution.  The  resulting  high  resistivity  layer  (at 
x  ~  lCji  in  Fig.  2C)  accounts,  of  course,  for  the  high  apparent  dc 
resistivity.  Capacitance  measurements  on  this  sample  would  yield  an 
anomalously  high  apparent  dielectric  constant,  approximately  in  the  ratio 
of  the  sample  thickness  to  the  thickness  of  the  insulating  layer. 

In  practice,  we  believe  that  this  situation  usually  prevails,  due 
to  heat  treatments  made  in  the  application  of  the  electrodes  or  other 
steps  in  the  preparation  of  the  specimen.  We  suspect  also  that  this 
phenomenon  is  not  confined  to  rutile  but  rather  may  occur  in  a  wide  class 
of  materials.  Note  that  the  only  requirements  are  a  work  function  greater 
than  that  of  the  electrodes  and  a  positive  ion  which  is  mobile  under  an 
electric  field  at  elevated  temperature,  but  not  at  room  temperature. 

In  general,  satisfactory  treatment  of  Eq.  (5)  requires  detailed1 
knowledge  of  trapping  levels  (which  have  been  ignored  at  present),  as 


well  as  precise  values  of  such  parameters  as  nQ*  These  calculations  are 
presented  solely  for  the  purpose  of  demonstrating  that  the  model  leads 
to  qualitatively  good  agreement  with  typical  experimental  results  to  be 
presented  in  the  next  section,  for  physically  reasonable  choices  of 
parameters. 

If  the  electrode  work  function  is  larger  than  the  rutile  work 
function,  electrons  will  flow  from  the  rutile  to  the  metal  until  the 
Fermi  level  is  constant  throughout.  This  will  result  in  an  (Ec  -  E^,) 
which  is  increased  at  the  surface  compared  to  the  value  before  electrodes 
were  attached.  In  the  middle,  (Ec  -  E^)  will  be  nearly  the  same  as  before, 
wit‘i  the  width  of  the  transition  region  being  determined  by  the  trap 
density.  If  the  sample  is  now  heated,  any  mobile  positive  ions  will  be 
drawn  toward  the  surface,  but  in  this  case  the  electric  fields  are  much 
smaller,  since  the  number  of  electrons  available  to  move  (<  nc  +  n^), 
and  hence  the  net  charge  density  producing  the  field,  is  much  smaller 
than  before.  We  have  carried  out  a  detailed  analysis  of  this  case, 
following  the  procedure  described  above;  unless  the  work  function  difference 
is  great  enough  to  induce  p-type  conductivity  (which  is  inconsistent  with 
experimental  results  except  for  certain  special  cases  to  be  discussed  in 
Section  III),  the  changes  produced  in  the  ion  distribution  do  not  result 
in  increased  conductivity  anywhere  in  the  crystal.  The  predicted  dis¬ 
sipation  and  capacitance  curves  are  the  same  as  those  in  Fig.  2A.  The 
corresponding  ion  distribution  is  a  slow  monotonic  decrease  from  the 
electrodes  toward  the  center  of  the  crystal;  the  resistivity  is  as  high 
or  higher  everywhere  as  it  was  before  application  of  electrodes.  The 
details  depend  on  trap  densities  and  initial  ion  concentration,  but  the 
general  result  is  as  outlined.  The  possibility  of  p-type  conductivity 


apparently  must  be  considered  in  special  cases,  but  cannot  usefully  be 
analyzed  on  the  basis  of  presently  available  data. 

The  third  case  (equal  work  functions)  does  not  involve  space 
charge  layers  and,  hence,  no  ion  motion  is  anticipated.  There  is  an 
additional  possibility,  that  the  electrode-rutile  interface  (or  the 
rutile  surface)  is  charged,  due  to  surface  states  or  defect-induced 
double  layers,  but  the  experimental  results  seem  to  be  adequately 
accounted  for  by  the  above  analysis. 

The  authors  are  not  aware  of  any  direct  measurements  of  the  rutile 
work  functionj crude  estimates  based  on  previous  calculations1^  of 
electrostatic  potential  of  charged  entities  in  the  rutile  lattice  indicate 
that  a  rather  large  work  function  should  be  expected.  This  is  particularly 
reasonable  in  view  of  the  large  polarizability  of  rutile.  If  the  assump¬ 
tions  of  our  model  are  correct,  the  existence  of  the  large  electronic 
dissipation  peaks  observed  in  this  system  would  indicate  a  work  function 
significantly  greater  than  the  ~  5  eV  which  is  typical  for  metals.  Our 

calculations  indicate  that  dielectric  loss  measurements  at  frequencies 

5 

<  10  Hz  are  insensitive  to  the  exact  value  of  the  work  function  of 
the  electrodes,  so  long  as  it  is  less  than  that  of  rutile  by  sw  0.15  eV. 


III.  Experimental  Results 

Measurements  of  C  and  D  have  been  made  on  a  variety  of  rutile 
crystals,  subjected  to  various  heat  treatments  and  surface  preparations, 
and  with  several  electrode  materials.  Measurements  reported  herein  were 
made  on  thin  disk-shaped  crystals,  oriented  with  the  c-axis  perpendicular 
to  the  face  of  the  disk  (ci).  A  number  of  measurements  were  also  made 
on  crystals  with  the  c-axis  parallel  to  the  major  face  of  the  disk; 
behavior  of  these  specimens  paralleled  that  of  the  Qt  specimens,  after 
allowing  for  the  smaller  dielectric  constant,  except  that  the  temperature- 
induced  changes  proceeded  more  slowly.  Typically,  equivalent  results 
could  be  obtained  in  comparable  times  by  raising  the  temperature  by 
~  100°C.  Sample  thickness  ranged  from  250  to  600|i.  Results  obtained 
were  not  significantly  dependent  on  sample  thickness,  although  longer 
heat  treatments  were  required  to  produce  the  same  results  in  thicker 
crystals.  All  measurements  were  made  using  a  General  Radio  Model  1615-A 
bridge.  Guard  electrodes  were  used  in  some  cases,  but  did  not  appear 
to  affect  the  measurements,  indicating  that  surface  currents  were  not 
Important .  "Pure"  crystals  from  various  sources  were  tested;  some 
differences  were  observed  but  there  did  not  appear  to  be  a  consistent 
correlation  with  source.  Doped  crystals  exhibit  quite  different  and 
variable  behavior,  but  not  enough  data  is  yet  available  to  justify 
discussion. 

The  results  obtained  on  various  sanples  and  electrode  materials 
were  widely  variable  and  obviously  dependent  on  a  large  number  of 
parameters,  not  all  of  which  can  be  adequately  controlled  or  evaluated 
with  presently  available  crystals.  For  this  reason,  it  is  not  possible 
to  present  a  conplete  and  unified  treatment  of  this  complex  problem  at 
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present,  and  no  attempt  will  be  made  to  present  data  in  detail;  rather, 
we  singly  describe  our  results  qualitatively,  drawing  conclusions  where 
possible,  and  indicate  areas  where  additional  study  is  required.  Our 
results  are  presented  in  terms  of  the  observed  effects  of  various 
experimental  parameters;  in  some  cases,  considerable  interpretation  is 
involved  and,  hence,  the  conclusions  reached  must  be  regarded  as 
tentative. 


A.  Surface  condition 

The  presence  of  surface  damage  results  in  complicated  behavior, 
which  has  not  been  carefully  studied.  In  general,  however,  the  results 
appear  to  be  consistent  with  the  assumption  that  surface  damage  produces 
a  layer  of  relatively  high  conductivity  in  which  diffusion  of  Ti 
interstitials  or  other  donor  centers  (discussed  below)  is  inhibited. 

We  were  unable  to  reduce  the  dielectric  loss  below  ~  0.1  until  after 
removal  of  mechanical  damage  by  etching,  or  after  annealing  above 
1000°C.  If  mechanical  damage  results  in  reducing  (Ec  -  E^,)  in  the 
damaged  region  (for  example,  by  introducing  interstitial  defects), 
these  observations  would  be  consistent  with  the  analysis  in  the  previous 
section,  with  the  loss  resulting  from  one  of  the  following  mechanisms: 

1)  In  the  absence  of  heat  treatments  the  loss  results  from  the  low 
resistivity  damaged  surface  layer. 

2)  If  the  crystal  is  heated  to  moderate  temperatures  the  loss  results 
from  the  induced  ion  migration  described  in  Sec.  II  with  the 
damaged  layer  playing  the  same  role  as  a  metal  electrode  (i-e.,  a 
material  with  higher  initial  Fermi  level). 

High  temperature  (~1000°C)  treatment  would  allow  annealing  of  these 
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defects,  resulting  in  a  raising  of  (E  -  E^)  in  the  damaged  layer  and 
reduction  of  the  induced  dissipation.  Further  evidence  of  the  sensitivity 
of  dielectric  loss  to  surface  damage  is  provided  by  the  observation  that 
mechanical  polishing  of  an  etched  crystal  with  O.ljj  alumina  produced 
easily  detectable  changes. 

All  results  described  below  were  obtained  on  crystals  which  had 
been  carefully  lapped,  followed  by  etching  in  KOH  at  ~  1+50°C  to  remove 
at  least  25  microns. 


B.  Electrodes 

A  variety  of  electrode  materials  were  applied  to  the  crystals, 
including  Ag  lacquer,  colloidal  graphite,  PT  and  Au  pastes  (fired  @ 
~900°C)  and  evaporated  Ti,  Al,  Ag,  and  Au.  All  electrode  materials 
which  could  be  applied  at  room  temperature  produced  indistinguishable 
results  on  a  given  sample  so  long  as  the  crystal  was  not  heated  after 
application  of  electrodes.  Loss  curves  for  different  sanqples  were  widely 
variable,  depending  on  heat  treatment  prior  to  application  of  electrodes, 
but  were  reproducible  for  a  given  crystal  with  different  electrodes. 

The  observed  dissipation  depended  in  a  complicated  way  on  sample  purity 
sis  well  as  temperature,  ambient,  and  cooling  rate  associated  with  previous 
heat  treatment.  The  mechanisms  responsible  for  the  loss  curves  observed 
are  not  well  understood;  the  only  thing  that  can  be  said  at  present  is 
that  the  loss  is  definitely  electronic  in  origin.  (See  Sec.  I.)  It  is 
worth  noting  that  severe  oxidation  (i.e.,  dry  Og  at  1000°C)  does  not 
always  produce  the  low  loss  which  would  be  expected  from  dc  conductivity 
measurements,  suggesting  that  it  is  probably  possible  (at  least  in  some 
crystals)  to  produce  weak  p-type  conductivity  at  room  temperature; 
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failure  to  observe  corresponding  dc  conductivity  probably  results  from 
an  intrinsic  layer  induced  by  the  presence  of  electrodes  (due  to  a 
transition  from  n-type  to  p-type  conductivity). 

Considerable  progress  has  been  made  in  understanding  the  nature 
of  the  changes  which  occur  when  crystals  are  heat-treated  with  electrodes 
in  place.  It  is  possible  to  produce  dielectric  loss  D  <  .01  over  the 
frequency  range  from  20  to  10  cycles  (with  constant  capacitance  equal 
to  the  calculated  value  to  within  1$,  assuming  e  =  173  eQ)>  or  loss 
peaks  anywhere  in  this  frequency  range  with  D  >  10. 

The  loss  curves  observed  in  crystals  heated  with  electrodes  in 
place  can  be  understood  in  terms  of  the  model  described  in  the  previous 
section.  However,  due  to  differences  in  sample  purity,  thermal  history, 
etc.,  the  data  are  so  widely  variable  (although  accurately  reproducible 
for  a  given  crystal,  electrode  material,  and  heat  treatment)  that  no 
useful  purpose  would  be  served  by  a  detailed  presentation  of  results. 
Figure  k  is  typical  of  crystals  exhibiting  moderately  high  loss,  in 
that  the  peak  shape  is  roughly  the  same  in  most  cases,  and  the  dissi¬ 
pation  drops  off  to  ~1Q$  of  the  peak  value  at  high  and  low  frequencies. 

The  electrode  materials  fall  into  two  categories;  those  which 
produce  irreversible  effects  on  heating  and  those  which  do  not.  By 
irreversible  we  mean  that  the  changes  in  the  loss  spectrum  produced  by 
a  given  heat  treatment  (at  temperatures  of  ^00°C  or  less)  were  not 
completely  reversed  by  removal  of  electrodes  and  reheating  to  the  same 
temperature.  All  of  the  electrode  materials  listed  above,  except  for 
Ag  and  Au,  resulted  in  irreversible  changes,  all  of  which  were  qualita¬ 
tively  similar  (although  the  temperature  dependence  was  different)  and 
could  be  accounted  for  by  diffusion  of  donors  into  the  crystal  from  the 
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electrode-crystal  interface.  These  materials  were  also  qualitatively 
different  from  Au  and  Ag  in  that  the  effect  produced  by  heating  did  not 
"saturate"  (by  "saturate"  we  mean  that  increasing  the  length  of  the  heat 
treatment  produces  no  further  changes)  while  changes  produced  by  Au  and 
Ag  did.  For  example,  a  500|x  cx  disk  with  Au  or  Ag  electrodes  saturates 
in  ~  10  hours  at  350°C,  while  irreversible  changes  associated  with  A1 
or  Ti  electrodes  are  not  observed  to  saturate  after  100  hours.  The 
irreversible  changes  apparently  are  produced  either  by  diffusion  of 
material  from  the  electrode  (Al,  for  example )  or  by  chemical  reaction 
of  the  electrode  material  with  the  rutile,  resulting  in  excess  Ti  which 
diffuses  into  the  crystal.  Al  appears  to  diffuse  into  the  crystal, 
presumably  as  a  3+  interstitial,  since  the  dissipation  increases  con¬ 
siderably  faster  in  this  case  than  with  Ti  electrodes.  C  electrodes 
were  not  investigated  in  detail,  except  to  establish  that  irreversible 
effects  occurred  at  400°C-  Pt  and  Au  paste  electrodes  require  firing 
at  elevated  temperatures  and  were  abandoned  early  in  the  investigation. 

The  behavior  of  crystals  with  evaporated  Ag  or  Au  electrodes 
apparently  is  quite  well  accounted  for  by  the  model  discussed  in  the 
preceding  section.  We  find  no  evidence  for  contamination  of  the  crystal 
or  reaction  with  the  electrodes,  and  all  changes  observed  are  fully 
reversible  and  reproducible  for  a  given  crystal,  although  the  observed 
dissipation  depends  strongly  on  the  temperature  and  the  ambient  atmos¬ 
phere  of  the  heat  treatment  (as  discussed  below),  as  well  as  on  high 
temperature  heat  treatment  prior  to  electrode  application  and  on  sample 
purity.  Thus  we  conclude  that  the  observed  dissipation  in  this  case 
results  from  the  nonuniform  distribution  of  donors,  induced  by  the 
difference  in  work  function  of  the  electrode  and  the  rutile,  as  described 
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in  the  preceding  section.  Au  and  Ag  electrodes  apparently  produce 
indistinguishable  results;  our  calculations  indicate  that  the  difference 
in  work  functions  for  Au  and  Ag  would  not  result  in  detectable  differ¬ 
ences  in  behavior  so  long  as  both  are  significantly  smaller  than  that 
for  rutile.  In  fact,  we  would  expect  virtually  identical  results  for 
any  nonreactive  metal  electrode. 

Results  obtained  with  Ag  lacquer  electrodes  are  similar,  though 
somewhat  more  variable,  but  again  are  fully  reversible.  It  appears 
likely  that  hydrogen  from  the  organic  binders  could  diffuse  into  the 
crystal  in  this  case;  our  observations  are  consistent  with  this  inter¬ 
pretation. 

A  common  denominator  in  all  the  above  results  is  the  critical 
dependence  of  room  temperature  dielectric  loss  on  any  heat  treatments 
given  the  sample  during  or  after  application  of  electrodes.  Heating 
to  100°C  for  a  few  minutes  frequently  resulted  in  changes  in  D  by 
factors  of  two  or  more,  and  minor  changes  were  even  observed  after 
"annealing"  at  room  temperature.  In  general,  large  loss  peaks  were 
observed  in  crystals  heated  to  250°C  or  more,  then  rapidly  cooled  to 
room  tenperature.  Prolonged  annealing  at  temperatures  below  250°C  or 
slow  cooling  (~10°C/hr  or  less)  resulted  in  low  loss,  although  results 
depended  strongly  on  ambient  during  heat  treatment,  sample  purity,  and 
thermal  history. 

Because  of  the  irreversible  changes  with  other  electrode  materials, 
the  remainder  of  our  discussion  will  be  confined  to  results  obtained 
with  Ag  and  Au  electrodes.  "Typical"  results  for  a  nominally  stoichio- 
metric  crystal  (p  >  10  Q  cm,  as  inferred  from  two  terminal  dc 

fV 

conductivity  measurements)  at  room  temperature  are  presented  in  Fig.  4, 
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in  terms  of  c/CQ  (where  C  is  the  measured  series  equivalent  capacitance 
and  CQ  is  the  calculated  capacitance  for  a  loss-free  crystal,  assuming 
c/eQ  =  173)-  The  crystal  was  500|i  thick,  had  been  carefully  lapped  and 
etched  to  remove  surface  damage,  and  annealed  in  air  at  800°C.  Ag 
electrodes  were  applied  at  room  temperature;  the  measured  dissipation 
and  capacitance  are  shown  as  Curves  A  in  Fig.  4.  Subsequent  heat 
treatment  (with  electrodes  in  place)  for  12  hours  in  air  at  350°C, 
followed  by  rapid  cooling  to  room  temperature  resulted  in  Curves  B. 

The  striking  changes  in  C  and  D  are  fully  reversible;  repeating  the 
350°C  heating  cycle  after  removal  of  electrodes  restored  the  behavior 
indicated  by  Curves  A.  It  was  not  always  possible  to  achieve  the  low 
loss  indicated  by  Curves  4-A;  however,  alternate  heat  treatment  with 
and  without  electrodes  always  produced  large  changes  in  C  and  D,  which 
were  reproducible  for  many  cycles.  Heating  to  lower  temperatures  with 
electrodes  and  quenching  results  in  a  shift  of  the  dissipation  peak 
to  lower  frequency  and  a  decrease  in  peak  height,  as  would  be  predicted 
by  the  analysis  in  Sec.  II. 

As  previously  noted,  the  observed  results  are  strongly  dependent 
on  a  number  of  experimental  parameters;  these  curves  should  only  be 
viewed  as  representative  of  results  obtained  for  the  conditions  specified. 

A  few  measurements  were  made  at  temperatures  intermediate  between 
room  temperature  and  annealing  temperature.  The  parameters  (C  and  D) 
generally  were  quite  unstable,  as  would  be  expected.  The  results  were 
in  qualitative  agreement  with  predictions  of  the  analysis  of  the 
preceding  section. 
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D.  Effect  of  Ambient  During  Heat  Treatment 
The  dissipation  curves  observed  for  a  given  heat  treatment  of  a 
particular  crystal  were  found  to  ,  depend  strongly  on  the  atmosphere  in 
which  the  treatment  was  carried  out.  The  results  obtained  appeared  to 
depend  only  on  the  equivalent  Hg  partial  pressure  of  the  atmosphere  and 
could  in  all  cases  be  interpreted  at  least  qualitatively  in  terms  of 
the  analysis  of  the  preceding  section,  by  assuming  that  the  migrating 
donors  were  H+  interstitials  and  that  the  equilibrium  H+  concentration 
in  a  given  crystal  was  determined  by  the  equivalent  Hg  pressure  in  the 
ambient.  Among  the  ambients  tested  (in  order  of  increasing  Hg  partial 
pressure)  were  Og  (with  liquid  nitrogen  cold  trap),  Og  (saturated  at 
room  temperature  with  HgO )  pure  HgO  vapor,  NH^-HgO,  and  pure  Hg 
(p  <  0.01  Torr).  The  general  behavior  for  crystals  quenched  from  a 
given  annealing  temperature  in  various  Hg  partial  pressures  was  for 
the  dissipation  peak  to  shift  to  higher  frequencies  with  increasing 
Hg  pressure,  as  would  be  expected  from  a  raising  of  the  Fermi  level 
in  the  crystal.  Thus,  atmospheres  with  very  small  partial  pressures 
of  Hg  generally  resulted  in  loss  peaks  at  low  frequency;  often  the 
peak  appeared  to  be  below  the  experimentally  accessible  range  and  only 
a  low  frequency  "tail"  was  observable. 

One  important  exception  to  this  behavior  was  observed.  A  250|u 
ci  crystal  of  relatively  high  purity  (no  quantitative  information 
available)  was  etched  and  strongly  oxidized  (Og,  100  Torr,  with  liquid 
nitrogen  cold  trap,  kO  hours  at  950°C,  slow  cooled)  before  applying 
Au  electrodes.  Subsequent  heat  treatment  in  dry  Og  at  35C°C  produced 
a  large  loss  peak  at  ~  5  kHz,  which  shifted  to  lower  frequency  and 
disappeared  at  successively  higher  IL  partial  pressure.  D  was  <  0.01 
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over  the  frequency  range  from  20  to  IQr  Hz  after  heating  in  20  Torr  of 
HgO.  Further  increase  in  Hg  pressure  resulted  in  the  usual  behavior 
(a  loss  peak  moving  to  higher  frequencies  as  Pu  increased).  This  is 
strongly  suggestive  of  p-type  behavior  (which  has  not  previously  been 
convincingly  demonstrated  in  TiOg),  produced  by  a  Fermi  level  near  the 
top  of  the  valence  band  in  the  bulk  of  the  crystel,  with  an  insulating 
layer  beneath  the  electrodes  resulting  from  n-type  electrode  contact. 

The  H  interstitial  concentration  in  these  crystals  was  independently 
monitored  by  observation  of  the  infrared  absorption  associated  with  the 
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OH  stretching  vibration.  Although  the  concentration  was  strongly 
sample-dependent  (presumably  due  to  variation  in  electron  trap  density), 
the  results  were  consistent  with  the  interpretation  above,  in  that  the 
H  concentration  increased  monotonically  with  increasing  Hg  partial 
pressure  of  the  annealing  atmosphere.  Measured  H  concentration  ranged 
fiom~  5x10  to  10  7cmJ  for  the  various  atmospheres  described  above. 

This  is  a  surprisingly  narrow  range  of  concentrations  in  view  of  the 
fact  that  (Hg  partial  pressure  in  the  ambient)  changed  by  many  orders 
of  magnitude  (the  equivalent  P  in  a  dry  0o  atmosphere  is  undoubtedly 

Hg  2 

determined  by  outgassing  of  hydrocarbons  from  o-rings,  etc.,  but  is 
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unlikely  to  be  larger  than  10  Torr).  However,  this  weak  dependence 

on  P„  is  readily  accounted  for  in  a  crystal  such  as  rutile  by  an 
"2 

analysis  of  the  thermodynamics  of  the  process.  It  is  readily  shown  that 


„  .  /  r,  -  E-/2KT 

cc  (P  )^  e  F 

where  n^+  is  concentration  of  H+  interstitials  and  is  the  electron 

Fermi  level  in  the  crystal.  E_  may  change  by  an  amount  comparable  to 

r 

the  band  gap  (~  3*2  eV  in  TiOg)  as  n^+  is  varied,  if  there  is  some 
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concentration  of  deep  electron  trapping  levels  in  the  crystals,  and  H+ 
is  the  dominant  donor  impurity.  This  readily  accounts  both  for  the 
limited  solubility  of  H+,  even  in  a  pure  Hg  atmosphere,  and  tne  extreme 
difficulty  in  reducing  n„+  below  ~  10  /cm  .  The  strong  dependence  of 

ri 

H+  solubility  on  E_  probably  also  accounts  for  the  difficulty  in  making 

r 

meaningful  measurements  of  H+  diffusion  in  TiOg. 

The  observed  dependence  of  P^.  appears  to  establish  the  F+  inter- 

stitial  as  the  required  mobile  impurity  in  at  least  some  cases.  Whether 

or  not  it  is  the  only  mobile  ion  is  a  question  which  cannot  be  answered 

on  the  basis  of  presently  available  results.  The  possibility  of  other 
/  44 

ions  (particularly  T1  ,  which  is  known  to  be  present  in  substantial 
concentrations)^  participating  in  the  process  cannot  be  ruled  out. 

The  equilibrium  configuration  for  two  or  more  mobile  ions  is  likely  to 
be  quite  complicated,  due  to  the  dependence  of  solubility  on  electro¬ 
static  potential  discussed  above.  The  rather  surprising  results  obtained 
from  calculations  for  the  single  ion  case  suggest  that  the  resistivity 
profile  resulting  from  redistribution  of  several  ion  species  could  be 
very  complex. 


E.  Voltage  Dependence  of  Measured  Dissipation 
In  most  cases,  the  measured  dissipation  was  essentially  independent 
of  the  amplitude  of  the  ac  voltage  applied  to  the  crystal,  at  least  up 
to  ~  30V  rms.  The  results  described  above  refer  only  to  such  cases,  or 
to  crystals  for  which  the  measured  dissipation  was  constant  for  sufficiently 
small  applied  voltage.  In  some  crystals,  however,  a  striking  voltage 
dependence  of  both  capacitance  and  dissipation  was  observed.  These  were 
crystals  which  had  been  equilibrated  in  rather  high  partial  pressures  of 
Hg*  presumably  resulting 
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in  a  Fermi  level  fairly  near  the  conduction  band.  In  some  cases,  the 

dissipation  at  low  frequencies  (<  10  Hz)  was  observed  to  change  by  as 

much  as  a  factor  of  10  as  the  applied  voltage  was  increased  in  the  range 

from  1  to  30  volts  rms.  A  similar  voltage  dependence  of  the  apparent  dc 

conductivity  was  also  observed.  This  effect  is  not  to  be  confused  with 
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the  time-dependent  conductivity  frequently  observed  in  rutile,  which 
is  presumably  due  to  field-induced  ion  motion.  Some  minor  changes  in 
our  observed  conductivity  occurred  several  hours  after  applying  a  dc 
field,  but  the  effect  described  here  was  clearly  different,  since  changes 
in  conductivity  occurred  on  a  millisecond  time  scale.  As  the  frequency 
was  increased,  the  voltage  dependence  of  D  diminished  and  was  essentially 
unobservable  at  10  Hz  at  voltages  up  to  100  volts.  The  observed  behavior 
is  well  represented  by  an  equivalent  circuit  consisting  of  an  ideal 
loss-free  capacitor  (  with  c  equal  to  the  calculated  capacitance  of  the 
crystal),  in  parallel  with  a  voltage-dependent  resistor.  Interestingly, 
the  frequency  below  which  D  increased  substantially  above  the  small 
signal  level  at  a  given  voltage  corresponded  approximately  to  the  recip¬ 
rocal  of  the  transit  time  for  an  electron,  predicted  from  the  applied 

2  3 

voltage  and  observed  carrier  mobility  (~  0.3  cm  /volt  sec).  This 
appears  to  imply  that  it  is  possib]e  for  carriers  to  be  injected  at  one 
electrode,  and  travel  completely  through  the  specimen  without  being 
trapped— which  would  be  surprising  to  say  the  least. 

IV.  Conclusions 

It  is  difficult  to  draw  many  firm  conclusions  from  the  admittedly 
qualitative  icsults  presented  here.  However,  the  results  are  sufficiently 
interesting  and  suggestive  to  merit  a  tentative  interpretation. 
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Part  of  the  observed  behavior--the  enhancement  of  dissipation 

which  results  from  heating  and  quenching  of  a  crystal  with  electrodes — 

appears  to  be  at  least  qualitatively  explained  by  the  electronic  loss 

mechanism  outlined  in  Sec.  II,  which  we  believe  has  not  been  previously 

suggested.  It  seems  likely  that  this  loss  mechanism  may  be  fairly 

common  in  insulating  crystals,  particularly  in  the  titanates,  and  in 

other  transition  metal  oxides.  Similar  though  less  dramatic  observa- 
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tions  have  been  reported  in  Al^O^  and  MgO,  neither  of  which  are 
particularly  good  candidates  for  such  a  mechanism,  because  of  their 
large  band  gap  and  small  e*  Other  observations,  such  as  the  dependence 
of  observed  loss  on  heat  treatment  prior  to  electrode  application, 
sample  purity,  and  surface  damage,  cannot  be  accounted  for  in  detail 
at  present,  but  are  probably  the  result  of  known  mechanisms. 

However,  some  of  our  observations  do  not  appear  to  be  consistent 
with  current  understanding  of  this  crystal,  and  merit  much  closer 
examination.  Of  particular  interest,  because  of  the  apparent  implica¬ 
tions,  are  l)  the  voltage  dependence  of  D  observed  in  some  cases,  and 
2)  the  large  loss  in  some  heavily  oxidized  crystals,  which  first 
decreases,  then  increases  again  after  heating  in  successively  higher 
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APPENDIX  A 


The  method  chosen  employs  the  Tchebycheff  technique  for  numerical 
integration.  This  is,  of  course,  not  the  only  possible  approach  and  in 
some  of  our  cases  not  the  best.  It  is  the  simplest  and  has  proven  adequate 
for  our  purposes. 

In  the  Tchebycheff  method  of  numerical  integration,  the  integral  is 
approximated  by:  . 


J  <Ku)du  =  “  r<J’(u1)  +  *(u2)  +  — +♦(%)] 


where  the  are  the  zeros  of  the  Tchebycheff  polynomials.  This  formula 


can  be  extended  to  an  arbitrary  range  of  integration  by  a  transformation  of 
variables : 


where: 


b  I 

V  y(x)dx  =  (b-a)^$(u)du 
a  ■+ 


x  =  (b-a)u  +  (^gS) 

♦(«)  =  y  f" (b-a)u  +  (■—)] 


Hence  the  expression  for  D(u))  becomes: 


n 

E 


Pb 


D(uu)  «  b=l  (p^oue)  +  1 


.  2 

n  eou 


b=l  (p-uite)  +  1 


We  then  form  the  function 

n 


♦  *  E  (log  D,  -  log  D  ) 

J=1  J  J 

A 

where  Dj  is  the  measured  dissipation  and  Dj  the  calculated  dissipation 
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at  ou  =  oUj*  We  then  choose  the  set  p^  to  minimize  A  computer  program 
has  been  written  to  accomplish  this.  It  calculates  the  required  change 
in  the  set  p^  assuming  the  original  set,  p^,  is  close  enough  to  the 
minimum  for  a  quadratic  fit  to  be  accurate.  It  then  calculates  the 
gradient  and  makes  a  step  whose  length  and  direction  is  a  compromise 
between  the  quadratic  and  gradient  methods.  The  relative  weight  given 
each  is  determined  by  the  angle  between  the  two  directions,  i.e.,  when 
the  angle  is  large  the  gradient  method  is  heavily  weighted.  Once  the 
optimum  set  pk  is  determined,  it  is  plotted  against  the  roots  of  the 
Tchebycheff  polynomials  to  obtain  the  resistivity  spectrum. 
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1. 


Figure  Captions 


a)  Assumed  band  configuration  of  isolated  rutile,  and  metal  electrode. 

b)  Band  configuration  of  system  after  electrical  contact  is  established. 


Location  of  conduction  band  minimum  relative  to  Fermi  level  vs  position 
in  crystal,  with  metal  electrode  at  x  =  0.  Curve  A  is  initial  config¬ 
uration  at  room  tenperature .  Curve  B  is  the  calculated  equilibrium 
configuration  at  350°C,  after  ion  migration.  Curve  C  is  the  corresponding 
configuration  after  quenching  to  room  temperature. 


3- 

Calculated  room  tenperature  capacitance  and  dissipation  for  band 
configurations  in  Fig.  2.  Curves  A  (this  figure)  corresponds  to 
Fig.  2A,  curves  B  correspond  to  Fig.  2C.  CQ  is  the  calculated  (loss- 
free)  capacitance  assuming  e  =  173  eQ* 

4. 

Measured  room  tenperature  capacitance  and  dissipation  for 
crystal  described  in  text.  Curves  A  were  obtained  before  heat 
treatment,  curves  B  after  equilibrating  at  350°C  in  air  and 
quenching  to  room  temperature. 
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Fig.  1 

a)  Assumed  band  configuration  of  isolated  rutile,  and  metal 
electrode. 

b)  Band  configuration  of  system  after  electrical  contact  is 
established. 
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C/C0 (B) 
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Calculated  room  temperature  capacitance  and  dissipation  for  band  configurations  in  Fig.  ' 
Cui-zes  A  ( this  figure )  correspond  to  Fig.  2A,  curves  B  correspond  to  Fig.  2C.  C  is  the 


